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The cyclotron resonance absorption of two-dimensional electrons in semiconductor heterostruc- 
tures in high magnetic fields is investigated. It is assumed that the ionized impurity potential is 
a dominant scattering mechanism, and the theory explicitly takes the Coulomb correlation effect 
into account through the Wigner phonons. The cyclotron resonance linewidth is in quantitative 
agreement with the experiment in the Wigner crystal regime at T = 4.2K. Similar to the cyclotron 
resonance theory of the charge density waves pinned by short-range impurities, the present results for 
the long-range scattering also show the doubling of the resonance peaks. However, unlike the case of 
the charge density waves, our theory gives the pinning mode independent of the bulk compressibility 
of the substrate materials. 
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I. INTRODUCTION 

Two-dimensional (2D) electron system in strong mag- 
netic fields has been studied intensively for many yearsi^. 
There are numerous experimetal and theoretical investi- 
gations concentrating on the low density electron sys- 
tem in high magnetic fields when the electron correlation 
becomes strong and the system can be considered as a 
Coulomb liquid or Wigner lattico?' 4 '^ 1 ?^!?!?. In order to 
study the characteristics of the 2D electrons in a mag- 
netic field, the cyclotron resonance (CR) is considered 
as one of the basic techniques, and important informa- 
tion about the 2D electron system can be obtained from 
the CR. Though there are many theoretical works con- 
sidering the CR in the highly correlated electron liquid 
and correlated electrons which form a Wigner crystal on 
the surface of liquid heliumi^S^, the CR of highly corre- 
lated electrons in semiconductor heterostructures is still 
an open question. 

In a magnetic field normal to the plane, the kinetic 
energy of 2D electrons is quenched and electrons exe- 
cute a cyclotron motion with orbit radius propotional 
to the magnetic length l c = y/hc/eB which decreases 
with increasing magnetic field B. Once B becomes suffi- 
ciently high so that l c is small compared with the typical 
mutual distance between electrons, the crystallization of 
the electron system into a Wigner crystal is expected if 
the electron density is sufficiently low. In semiconduc- 
tor heterostructures, when the mutual Coulomb energies 
are dominant, a Wigner crystal or a charge density wave 
may be formed in the quantizing magnetic field limit. 
Recently, an insulating state was observed on monolayer 



systems 



such a state is interpreted as a Wigner 



crystal pinned by the disorder in the sample^. 

The early theory by Fukuyama and Lesi&, which con- 
sidered the pinning of the charge density waves of the 2D 
electrons in heterostructures under the influence of the 
magnetic field, has predicted that the absorption line- 
shape of the system has double resonance peaks: one is 
located at the cyclotron frequency and the other related 
to the impurity pinning. In this theory, the magnetic 



field is treated classically and the impurity potential is 
assumed as short-ranged ones, and the pinning mode de- 
pends on the bulk compressibility of the substrate mate- 
rials. 

In this work, we present the theoretical analysis of the 
dynamic structure factor (DSF) and the linewidth of the 
CR due to ionized impurity scattering of the 2D electron 
systems in quantizing magnetic fields. Electrons are as- 
sumed to form a Wigner crystal in which the electrons 
are localized and oscillate around their equilibrium posi- 
tions, and the electron correlation is taken into account 
through the Wigner phonons. The numerical results of 
the CR linewidth are compared with experiments in the 
ALrGai^As heterostructures 

The paper is organized as follows. In Sec. [HJ the ba- 
sic relations for the conductivity, linewidth and DSF are 
presented. In Sec. IIIII the expression and the results of 
numerical calculation of the DSF are given. We then 
evaluate the CR linewidth in Sec. IIVI A comparision 
with the experiment for density denpendence of the CR 
linewidth is also given in this section. The double peak 
structure similar to the case of the pinning of charge 
density waves^ of the diagonal conductivity is shown 
in Sec. [3 Finally, we summarize our results and give 
concluding remarks in Sec. IVII 



II. BASIC RELATIONS 

The absorption lineshape of the 2D electron system 
in a magnetic field is characterized by the conductiv- 
ity <j±(lj) — u xx (lo) ± ia X y[ui) which may be expressed 
through the memory function M(uj) by: 
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(1) 



where uj c = eB /m*c is the cyclotron frequency, — e, to*, n e 
are the charge, effective mass and 2D density of the elec- 
trons, respectively. 

We assume that the 2D electrons forms a Wigner crys- 
tal and the 2D position vector of the nth-electron is given 
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by r„ = R n + u„ , where R„ is the lattice site vector and 
u„ is the deviation from R„. Further the scattering of 
electrons is caused by the static impurity potential V(r) 
which is long-ranged: 



III. DYNAMIC STRUCTURE FACTOR 



A. Wigner phonon spectra 



x e 

V ^ = ^ K [(r-R. t ) 2 + z2]i/2' 



(2) 



{Ri} 



where (R.;, zi) is the coordinate of the zth-impurity which 
is distributed randomly, and k the dielectric constant of 
the substrate matter. In this case, the density-density 
correlation function S(q,t) is defined as 

S(q, t) = e lq ' R < e - l i- u ™( t )e- iq - u - (0) >, (3) 

R 

where R = R„ — R m and < ■ ■ • > denotes the thermal 
average. 

The memory function M{oj) is related to the density- 
density correlation function by 
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where [3 is the inverse temperature, A the area of the 
electron system, and S{(\, oj) is the Fourier transform of 
S(q,t) 



S{q,u) 



+ 00 
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and will be called the dynamic structure factor (DSF) 
hereafter. In eq. (@J, U q is the Fourier transform of the 
random averaged correlation function of ionized impurity 
potentials U(v) = < V(r)V(0) > r given bj*ii 



r,= (^ 2 ^erfc(d/A). 



(0) 



Here, rii denotes the impurity concentration and the im- 
purity distribution is characterized by < Zi > r = —d and 
< (zi +d) 2 > r = A 2 /2. 

The imaginary part of the memory function M(u>) will 
be called the width function and denoted by j(u>) since 
it characteries the absorption line shape. We have: 

7H = ImAf(c) = -Y^U^ReSiq^). 

q 

(7) 

In order to evaluate this quantity, the calculation of 
the DSF (J5J will be given in the next section. 



In the absence of the magnetic field and impurities, 
Wigner phonon frequencies in the long-wavelength limit 
are approximately given byi& 



(0) 
kZ 



(longitudinal) 



,(0) 



y/2 k 

ujq — , (transverse) 
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where uj — 2(7m e ) 3 / 4 (e 2 / kto*) 1 / 2 and k a — ^/47rn e . 

In the presence of ionized impurities, the Wigner 
phonon spectra may have a gapii 

,(°) 2 i „,2 



-kA 



where gap v is given by 
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in the longwavelength limit. Here e A (k) is the polariza- 
tion vector of the Wigner phonons with wave vector k 
and mode A in the absence of magnetic fields: 

ej(k) = p £t = e,x £| (k), (12) 

where e z is the unit vector normal to the electron plane. 

When a magnetic field is applied normal to the electron 
plane, it is well-known that the longitudinal and trans- 
verse modes of the Wigner phonons are coupled and split 
into two modified modes with frequencies f^kjC? = ±) 
given by 
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k± 



\/tJ 2 + (Uki + LU kt ) 2 ± ^UJ 2 + (w k ( - ^kt) 2 



(13) 

The Fourier transform of the displacement vector u n 
and the real momentum n„ are given, respectively, by 




E^k^+E^k)^. 
e z x u(k) 



, (14) 



E^k^-E^k^ 1,(15) 



where p(k) is the Fourier transform of the canonical mo- 
mentum, &kj and 6^ are the annihilation and creation 
operators of the phonon, respectively, and 



E±(k) 



e ; (k)=F« 



^k+ ^k 



^k+ ^k- 



e t (k). 



(16) 
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B. Dynamic structure factor 



In order to obtain the expression of the DSF, we cal- 
culate the density-density correlation function defined by 
eq. ©. We assume here that the correlation between 
two electrons at the same site is much larger than that 
of electrons at different sites and disregard the contribu- 
tions with R ^ in S(q,t) (single-site approximation). 
Using eq. lfH|) we get 



where 



5(q,*)w< e - iq - u »«e iq - u ™ (0) > 
= exp [W(q,t)-W(q,0)}, 



j=± 



(17) 



(18) 
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E 



q-E,(k) 



^[{l+n^e-^ + n^e^*] , (19) 

where N is the number of the lattice sites and rikj the 
Bose-Einstein distribution function of the (kj)-phonon. 

When the magnetic field is very high such that 
(3kuj c 3> 1 and uj c » cj , we have f2£ + ~ u)% + v 2 » w 2 

2 k- ~ 



and f2?, w wj^wL/u; 2 <C /3 2 , and therefore 



W + {q,0)-W + (d,t)=x(l-< 
xT 2 



—iu) c t 



), 



(q, 0) - W_ (q, t) = ^ (t 2 - ih/ti), 
where x = hq 2 /2m*u) c and T is defined by 
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with Ti being the collision time due to the impurity scat- 
tering. 

After performing integration with respect to k and t, 
we obtain the final expression of the DSF— 
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It is interesting to see that the structure of the DSF 
(|26[1 is quite similar to the DSF of the highly correlated 
electron liquid by Monarkha et a/. 9 . However, it differs 
from ours in that the expression of r c in ref. 9 corre- 



^2 
c — c ■ 



sponds to our expressions by r c = x y r 2 _ ; + xT 2 

Monarkha et al. introduced the correlation effect phe- 
nomenologically by assuming the Doppler shift due 
to the ultra-fast electrons. Moreover, the factor 
exp [f3h(uj — nuj c )/2 — x(Yf3/A) 2 ] is missing in their re- 
sults, leading to an overestimation of the contributions 
from the higher Landau levels n > 1. 
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and will be called the broadening parameter of the DSF. 

From eq. (|1U|) . T 2 can be expressed as the sum of two 
contributions from the electron-electron interaction and 
from the electron-impurity scattering: 



(23) 



FIG. 1: Contour map of S(q,ui) using material con- 
stants of GaAs/AlGaAs at B = 7.4 T, T = 4.2 K and 
n e = 2 • 10 10 cm -2 . Impurity concentration is set to be 
m = O.lnii, and d/A = 1.70. 



where r r _ e is related to the electron-electron interaction: 
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and r o _i describes the single-electron effect related to the 
electron-impurity scattering: 
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In the limit q — > 0, the DSF (|26|l has a singularity 
S(q, us) oc 8{uj) which is consistent with Kohn theorem 20 . 
At q 7^ 0, the DSF has a series of peaks at points 
[u> « nuj c + xT 2 f3/4:h, x « n— 1/2) as shown in Fig.^and 
the peak widths are broadened not only by the electron 
interaction but also by the impurity scattering through 
the parameter T. The peak height deceases with increas- 
ing x and becomes very small at i > 1. Hence the most 
importance contributions into DSF come from the terms 
with n < 1. 
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In the extremely high magnetic fields and lu = uj c 
(cyclotron resonance regime) when T can be considered 
very small, all electrons at the lowest Landau level, and 
thus the term with n = 1 dominates the DSF. In this case, 
our theory reproducts the results obtained previously by 
Dykman 3 . Otherwise, the contributions of higher Landau 
levels must be taken into account. 



decreases with increasing electron density. Our the- 
oretical result is in quantitative agreement with the 
experiment in the region of the electron density 
2 • 10 10 cm" 2 < n e < 6 • 10 10 cm" 2 (the filling factor 
v = 2Trn e (cH/eB) varies from 0.08 to 0.34), where the 
electron system is expected to crystalize into a Wigner 
crystal. 



IV. CYCLOTRON RESONANCE LINEWIDTH 



Using eq. I|26|l . we rewrite the CR linewidth of the 2D 
correlated electrons under quantizing magnetic fields in 
more explicit form 



r2_i . pro* 

7(w c ) = . _ smn 



HT 2 
xY- -F n 



{2£?uj c | 1 - 



, (27) 



where £ = y/l + (T/3/4) 2 and the function F n (z) is de- 
fined by 



F n {z) 



_ J & 



(for n = 1), 



-=z n+1 '*K n+i (z), (forn^l). 



(28) 
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FIG. 3: 7(oj c )/cj c vs. n e at B — 10 T and for several tern- 
Here, K v {z) stands for the modified Bessel function of Peratures T = 500 mK (solid curve), 150 mK (dashed curve), 



the second kind. 



and 100 mK (dotted curve) 



80 




1 1 1 1 1 1 1 

2 3 4 5 6 7 

n e (10 10 cm" 2 ) 



FIG. 2: 7(<^c) vs. electron density for the 2D electrons in 
heterostructure GaAs/AlGaAs at T = 4.2 K, B = 7A T. The 
crosshairs are taken from the experimental results of Chou et 
al&. 



Figure [3] illustrates the density dependence of the CR 
linewidth at lower temperatures. It is easy to see that 
the density denpcndence of the CR linewidth is sensi- 
tive to the change of the temperature. At relatively 
high temperatures (T > 500 mK), the CR linewidth de- 
creases monotonically with n e . For lower temperatures 
(T < 150 mK), the CR linewidth has a shallow minimum 
around n e w 3- 10 10 cm -2 and then increases slightly. At 
sufficiently low temperatures (T < 150 mK), where r e _ e 
can be considered smaller than r e _; due to impurities, 
so that r ps r c _j, the role of impurities becomes most 
important. 

Figure 21 shows the magnetic field dependence of the 
ratio ^(lo c )/lo c at n e = 3 • 10 10 cm~ 2 for different tem- 
peratures. The different behavours in the magnetic field 
dependence of the CR linewidth at several values of tem- 
peratures are clearly seen in this figure. At sufficiently 
high temperatures (T > 200 mK), -f(uj c )/uj c increases 
monotonically with B. At T = 200 mK, j(oj c )/oj c 
first decreases, and then slightly increases with B. At 
T = 100 mK, j(oj c )/u! c becomes nearly constant for 
large B (B > 9 T). This means that CR linewidth is 
propotional to B at T w 100 mK. For T < 100 mK, 
7(w c )/w c decreases for a whole range of B. 



The electron density dependence of the CR linewidth 
is shown in Fig. [21 The CR linewidth monotonously 
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FIG. 4: y(Uc)/u)c vs. magnetic field at n e = 3 ■ 10 10 cm -2 for 
several temperatures: 500 mK (dotted cuve), 200 mK (dashed 
curver), 100 mK (solid curve) and 50 mK (dash-dotted curve). 



filling factors v. Two peaks are visible in the conductiv- 
ity. It is well-known that Wigner phonon spectra in a 
magnetic field split into two modified modes f2± as given 
by eq. (13). In a high magnetic field case, from eq. (13) 
we have J7 + rs uj c and f2_ rs v 2 /u> c w (27T7i) . The peak 
at lo = uj c corresponds to the ordinary CR, and one at 
lower frequency lu — f2_ is related to the pinning mode 
due to impurities. This behavour is similar to the case 
of the charge density wave pinned due to impurities in 
the presence of magnetic fields. Fukuyama and Lee de- 
veloped the CR theory for the charge density waves with 
the short-ranged pinning potential^. The magnetic field 
was treated semiclassically and their result for f2_ is de- 
pends on the bulk compressibility of the material. In our 
case, the pinning is determined by the ionized impuri- 
ties which is long-ranged, and the pinning mode 0_ is 
independent of material constants of the semiconductor 
substrate where electrons are located. 



VI. CONCLUSION 



V. DOUBLE PEAK STRUCTURE OF THE 
CONDUCTIVITY 

Next, we turn to the diagonal conductivity. It is de- 
termined from eq. (1) by 

ff(w) = i(o+(w)+o-_(w)). (29) 
Figure [S] exhibits the real part of the diagonal conduc- 
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FIG. 5: Frequency dependence of the diagonal conductivity 
for several filling factors at B = 6.15 T and T = 1 K. 



tivity Hea xx (u>) as a function of lo/lo c for several Landau 



We have presented in this paper a theoretical analy- 
sis of the cyclotron resonance of the 2D electron system 
in heterostructures under quantizing magnetic fields at 
finite temperature. We studied the DSF and the CR 
linewidth when the electron correlation effect is taken 
into account through Wigner phonons. We have ex- 
pressed the DSF as a summation over all Landau lev- 
els and obtained the analystic expression of the CR 
linewidth. By taking the electron correlation into ac- 
count through Wigner phonons, the single-electron and 
many-electron effects are considered simultaneously not 
phenomenologically as was done in previous theories^. 
The electron density dependence of the CR linewidth is 
shown to be in the quantitative agreement with exper- 
imental results by Chou et al& in the regime where the 
crystallization of electrons is expected. 

Similar to the case of the charge density waves with the 
short-ranged pinning potential^, our results of the con- 
ductivity also show the doubling of the resonance peaks. 
But unlike the charge density wave case, the pinning 
mode in our theory is independent of the bulk compress- 
ibility of the substrate materials. 

We have studied also other range of parameters than 
experiments. The numerical calculations shows that the 
CR linewidth monotonously decreases with increasing 
electron density at sufficiently high temperatures. On the 
other hand, the density dependence at low temperatures 
(T < 150 mK) is quite different: a sharp decrease in the 
very low density region (n e < 2 • 10 10 cm -2 ) and a slight 
increase for higher densities. The dependence of the CR 
linewidth on the magnetic field is also evaluated. Specif- 
ically, at sufficiently high temperatures (T > 100 mK), 
the cyclotron resonance is broadened considerably when 
the magnetic field becomes high, and at lower tempera- 
tures the broadening of the CR peak becomes very small. 
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